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We study the collective motion of atoms confined in an optical lattice operating inside a high 
finesse ring cavity. A simplified theoretical model for the dynamics of the system is developed 
upon the assumption of adiabaticity of the atomic motion. We show that in a regime where the 
light shift per photon times the number of atoms exceeds the line width of the cavity resonance, 
the otherwise tiny retro-action of the atoms upon the light field becomes a significant feature of 
the system, giving rise to dispersive optical bistability of the intra-cavity field. A solution of the 
complete set of classical equations of motion confirms these finding, however additional non-adiabatic 
phenomena are predicted, as for example self-induced radial breathing oscillations. We compare 
these results with experiments involving laser-cooled 85 Rb atoms trapped in an optical lattice inside 
a ring cavity with a finesse of 1.8 x 10 5 . Temperature measurements conducted for moderate values 
of the atom-cavity interaction demonstrate that intensity-noise induced heating is kept at a very 
low level, a prerequisite for our further experiments. When we operate at large values of the atom- 
cavity interaction we observe bistability and breathing oscillations in excellent agreement with our 
theoretical predictions. 

PACS numbers: 32.80.Pj, 42.50.Vk, 42.62.Fi, 42.50.-p 



I. INTRODUCTION 

Atoms regularly spaced in optical lattices are a widely 
studied model system of quantum optics Q which re- 
ceives growing attention also from other areas of physics 
ranging from the fields of condensed matter to quan- 
tum information processing 3j . Magneto-optic trapping 
techniques typically provide optical lattices with an av- 
erage inter-atomic distance of several microns, where 
dipole-dipole or spin interactions are negligible. While 
in this case the atomic dynamics is well understood it 
appears to be of limited relevance for most conceivable 
applications. Schemes that provide mutual interactions 
among the atoms should significantly increase the use- 
fulness of optical lattices. Recently, the formation of 
optical lattices inside optical resonators has become a 
subject of extensive research 0,|HIS0> because cavity- 
mediated long range interactions arise, possibly useful 
for new quantum c omp uting schemes [8( or novel laser 
cooling methods 0, EI EB [3 which do not rely on 
cyclic spontaneous emission and thus apply to a wider 
class of species without degradation to be expected at 
high densities. 

The key to exploiting these new concepts is a profound 
understanding of the interaction of trapped particles with 
intra-cavity light fields. For single atoms interacting with 
only a few photons in cavities with a small mode volume 
well below 10 mm 3 studies have been successfully per- 
formed within the last few years 0,^3- However, in 
these experiments only a few atoms can be trapped for 
no more than a few 10 ms due to heating caused by inten- 
sity fluctuations of the intra-cavity field. More recently 
it has been pointed out, that high finesse cavities with 
large mode volumes exceeding 1 mm 3 may be utilized to 
confine large atomic ensembles with the advantage of far 
longer trapping times. In this case the high frequency 



noise components of the intra-cavity light field are sig- 
nificantly suppressed owing to the long lifetime of the 
photons £j . Although operation at large detunings from 
atomic resonance is required in order to prevent unde- 
sirable spontaneous photons, for sufficiently high values 
of the cavity finesse the strong coupling regime can be 
reached. In this case the collective interaction strength 
given by the light shift per photon times the number of 
trapped particles exceeds the resonance line width of the 
cavity, and thus the otherwise tiny retro-action of the 
atoms upon the light field becomes a significant feature 
of the system. 

In this paper we experimentally and theoretically ex- 
plore the motion of atoms trapped in an optical lat- 
tice formed inside a high finesse ring-cavity with a large 
mode volume. We discuss two different regimes of opera- 
tion characterized by the strength of the collective atom- 
cavity coupling. Our observations for moderate coupling 
strengths show that the experimental challenges intro- 
duced by the small cavity resonance bandwidth can be 
handled and a stable lattice can be produced. At large 
atom-cavity couplings where the backaction of the atoms 
upon the lattice potential becomes significant, we find 
non-linear dynamics of the intra-cavity field and a col- 
lective character of the atomic motion. In this regime 
phenomena like dispersive optical bistability and self- 
induced breathing oscillations arise. Non-linear dynam- 
ics due to optical pumping and saturation has been ob- 
served in low-finesse resonators operating close to an 
atomic resonance |lGl . Only recently absorptive optical 
bistabi lity was observed in a small-volume high-finesse 
cavity |17| . Our lattice operates far from resonance where 
optical pumping and saturation are negligible and the 
atoms merely act as a dispersive medium. 

The paper is organized as follows. In Sec. [n]we present 
a simple physical picture of the atom-cavity coupling in 
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terms of coherent forward and backward Rayleigh scat- 
tering. In Sec. 11111 we briefly discuss the complete set of 
semiclassical equations of motion of the system. These 
equations are difficult to solve leading us to develop a 
simplified but surprisingly accurate model of the com- 
plex system dynamics based upon the assumption that 
the atomic sample adiabatically adjusts to the poten- 
tial position and depth. A steady-state analysis for the 
phase and amplitude of the intra-cavity field predicts 
dispersive bistability. In Sec. II VI we discuss the charac- 
teristics of our experimental setup. Both cavity modes 
are externally fed by adjustable amounts of power from a 
single-mode laser diode. A fast servo-lock keeps the laser 
frequency in resonance with one of the travelling wave 
modes of the cavity. The properties of this lock and the 
implications for parametric heating and expected trap- 
ping times are discussed in Sec. [V] Noise measurements 
of the light transmitted through the cavity together with 
lifetime and temperature measurements of the trapped 
atomic ensemble show that undesired intra-cavity inten- 
sity noise can be maintained at extremely low values 
despite of the small cavity bandwidth. In Sec. IVII we 
present experiments in the strong coupling regime. For 
asymmetric pumping of the cavity the unlocked intra- 
cavity field exhibits distinctive non-linear behavior re- 
sulting in bistability of the intra-cavity field. We find 
an excellent agreement with the predictions of the adi- 
abatic model of Sec. IIIII In Sec. IVIII we explore non- 
adiabatic aspects of the atomic motion. Observations of 
radial breathing oscillations are discussed and a numeri- 
cal simulation based on solving the full set of equations 
of motion is presented. 



II. MODE SPLITTING 

The interaction of the atomic sample with the intra- 
cavity light field can be understood in terms of coherent 
Rayleigh scattering. Consider the two degenerate, mutu- 
ally counterpropagating, travelling wave modes of a ring 
resonator (discriminated by indices (±) in the following) 
with a common resonance frequency lo c , as depicted in 
Fig. ^i). Let us add atoms into the common mode vol- 
ume with a resonance frequency uj a >> uj c such that 
their interaction with the intra-cavity fields has entirely 
dispersive character. If the atomic sample is homoge- 
neously distributed it will merely give rise to Rayleigh 
scattering in the forward directions described by a con- 
stant index of refraction n± = 1 + N Aq/lu Ci where Ao 
is the light shift per photon and N is the number of 
atoms. This causes a common frequency shift N Aq for 
both modes [Fig. EJj)]. Assume that both modes are ex- 
ternally coupled and an optical lattice is formed inside 
the cavity which tightly traps the atoms in the Lamb- 
Dicke regime, where the atomic distribution in each lat- 
tice site is confined to a fraction of the optical wavelength. 
In this case additional backscattering arises, which cou- 
ples the two counterpropagating modes and lifts their 




FIG. 1: a) Sketch of the experimental setup. Adjustable laser 
powers xo± are coupled into the two counterpropagating trav- 
elling wave modes of a ring resonator. The power transmitted 
through the cavity can be detected by a photo diode (PD). 
b) Mode splitting mechanism (explained in the text). 



frequency degeneracy giving rise to a frequency split- 
ting 2A^Ao \g\. The degree of atomic localization is mea- 
sured by the parameter g = -L e~ l2kZv ~ 2 ( x » +Vv ^ w ° , 
where k is the wave number, x v , y v and z v denote the 
atomic position coordinates, and wq is the mode ra- 
dius. If there is no statistical correlation between the 
axial and radial coordinates and the spatial distributions 
near each lattice site follow a Gaussian, we may write 

g= ^1 + 4 (oy/wo) 2 ^ exp(— 2k 2 a 2 ) exp(~ 2ikz cm ), i.e., 
the complex phase scales with the axial center of mass 
coordinate z cm = -k ~Y^ =1 z v and the modulus is a 
measure for the axial (cr z ) and radial (oy) spread of 
the atomic sample. Perfect localization corresponds to 
\g\ = 1, whereas a homogeneous distribution is described 
by \g\ = 0. 

If the interaction strength N A exceeds the cavity 
resonance bandwidth, the eigenmodes of the cavity ac- 
quire orthogonal standing wave geometries. One of them 
supports the optical lattice with atoms localized at the 
antinodes, i.e, the interaction is maximized and the corre- 
sponding frequency shift oj c + NAq (1 + \g\) exceeds that 
found for the travelling wave eigenmodes in presence of 
a homogeneous atomic sample. The second eigenmodc 
is not externally coupled and its nodes coincide with the 
atomic center of mass positions thus minimizing interac- 
tion. This leads to a reduced shift uj c + NA (1 — \g\) 
(Fignj>). The existence of a second empty eigenmode, 
blue-detuned with respect to the optical lattice, can be 
utilized for implementation of a sideband cooling scheme 
as described in Ref. [T^ . 

It is instructive to consider the change of the indices of 
refraction n± experienced by the travelling wave modes, 
due to back-scattering. According to the considerations 
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ofRef. [18f (Eg. 7) , to first scattering order, the modified 
refractive indices are given by n± = 1 +iV ^ (1 + 1 <7| | 2i | ) . 

Note that n + and n_ may differ, if the field amplitudes _ | y 

of the travelling waves a± are not equal, e.g., for asym- "' I 
metric external pumping of the cavity. This gives rise to 
interesting non-linear dynamics, if back-scattering be- 
comes relevant. A small external pumping asymmetry 
can yield a relative change of the refractive indices and a 
corresponding spatial phase shift of the optical standing 
wave potential. As the laser frequency is kept actively in 
resonance with say the (+)-wave, the (— )-wave is tuned 
slightly out of resonance and the potential well depth 
decreases. This in turn decreases the degree of localiza- 
tion \g\ which reduces the phase shift. As we will see in 
the following, in a certain parameter range this dynamics 
provides more than one steady-state, thus giving rise to 
bistability phenomena. 
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III. ADIABATIC MODEL 

We restrict ourselves to the case where the cavity 
modes support coherent fields with large mean photon 
numbers interacting with large thermal atomic samples 
at sufficiently high temperatures such that a classical de- 
scription of the optical intra-cavity fields and the atomic 
position and momentum variables is appropriate. More- 
over, we assume a large detuning of the pump frequency 
from the atomic resonance and hence neg ligible satura- 
tion. In this case according to Ref. [ijj the following 
equations describe the time evolution of the complex 
fields a+ and a_ for the counterpropagating travelling 
wave modes scaled to the field per photon 



dt 
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Here r) + and r/„ are the complex amplitudes of the in- 
coupled light fields, and 5 C denotes the detuning of the 
incoupled frequency from the resonance frequency uj c of 
the empty cavity. The radial bunching parameter g r ac- 
counts for the radial position dependence of the atom- 
light coupling. The diagonal elements of the matrix M 
comprise the forward scattering terms NA$g r which only 
depend on the degree of radial bunching while the off- 
diagonal terms which act to couple a+ and a_ are due 
to backward scattering and thus involve the localization 
parameter g. 

The equations for the momenta p v = (f> x ,„, Vy,v, Pz,u) 
and the positions x v = (x u , y u , z„) read 



(2) 

In total we are concerned with 6N + 2 coupled first- 
order non-linear differential equations for the 3N posi- 
tion and momentum variables, respectively and the two 
complex field amplitudes for the counter-propagating 
cavity modes. 

In experiments with small bandwidth cavities the fre- 
quency of the incoupled laser beam needs to by actively 
controlled in order to maintain some resonance condition, 
which allows to couple light into the cavity at all. In our 
implementation we servo lock the laser frequency to one 
of the counterpropagating travelling waves (say the (+)- 
mode) with a technique discussed in detail in [VJ Here 
we merely account for the fact that this servo acts to 
maintain the complex field amplitude a+ at a constant 
value, i.e., = 0. Hence, we may combine the two 

field equations JIJ into a single equation for the unlocked 
mode a_ 



d .A N , 

— a- = i ga_ 

dt e + 



a_7c - iA Q Ng*e + + S-j c . (3) 



using the abbreviations e± = ^r]± for the steady state 
intra-cavity fields in absence of atoms. We may choose 
e± to be real without loss of generality, thus fixing the 
spatial phase of the intra-cavity standing wave for the 
case when no atoms are present. 

For numerical processing and comparison with ex- 
periments it is useful to work with unit-free quanti- 
ties. We hence define scaled fields a = a_/-\//o and 
y/Xo± = £±/V^o with lo = |e+| 2 + |e-| 2 being the sum of 
the steady state intensities of each travelling wave mode 
in absence of atoms. Using the scaled time r = 7 c t and 
the scaled interaction strength U = Aq/j c we obtain 



dr 



UN 

Vxo+ 



\g\ga 2 -a-iUN\g\g*^/x^:+, 



(4) 



In order to incorporate the equations of motion for the 
atomic ensemble into this equation we seek expressions 
for the complex phase g = g/\g\ of the localization pa- 
rameter g, which is connected to the atomic center-of- 
mass coordinate, and the modulus \g\ which represents 
the spatial spread of the atomic sample in each poten- 
tial well. Therefore, let us consider atoms tightly con- 
fined in the optical lattice, such that the time scale of 
the axial motion is much shorter than the photon life 
time (l/27c) _1 which represents a lower bound for the 
time scale relevant for changes of the intra-cavity fields. 
In this case it is reasonable to expect, that the center- 
of-mass in each potential well of the optical standing 
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wave adiabatically follows the potential minimum, which 
formally is expressed by the relation g* — a/\a\. More 
specifically, we assume that the atomic distribution in the 
lattice can be described by a sum of Gaussians which are 
centered at the instantaneous antinodes of the lattice. To 
evaluate \g\ we take into account that there is no statis- 
tical correlation between the axial and radial coordinates 
writing \g\ w (e~ l2kz )(e~ 2r ^ w °), where the brackets de- 
note a Gaussian average. These expressions are readily 
calculated to be ( e - l2kz ) = e - 2k ^l and (fT 2 ^ K) = 
1/(1 + 4of /w 2 ). In order to determine the time evolu- 
tion of ay and a z we assume that the thermal atomic 
sample adjusts adiabatically to the potential well depth 
keeping the Boltzmann factors ^ ax = kBT ax /u> Viax and 
^rad = kBT ra d/LL> v ,rad constant, where ks is the Boltz- 
mann constant and T ax and T ra d are the axial and radial 
temperatures of the sample and U! v<ax and uj Vira d denote 
the axial and radial trap frequencies. By harmonically 
approximating the potential in axial and radial direc- 
tions we find 2a 2 /k 2 = k B T(t)/V Q . ax = t] ax y/a(0)/a(t) 
and 4a 2 /w 2 = r] ra d(^/XQ+ + |a(0)|)/(^/xoT + |a(f)|), with 
i] ax and r] ra d being the ratio between the thermal energy 
and the potential depth at t=0 for the axial and radial di- 
rections, respectively. With these approximations we are 
in the position to state a single equation for the complex 
unstabilized electric field amplitude: 



dr 



■UN _. .. 

i M a ) l a l a — a + 
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In this equation all the parameters can be easily ob- 
tained by measurements and a numerical integration 
yields simulations of the intra-cavity field, that can be 
compared to experimental data. In Sec. I VII we will show 
that despite its simplicity the adiabatic model presented 
here reproduces our observations very accurately. 

We may obtain additional physical insight into the dy- 
namical properties of Eq. JSJ) by representing the com- 
plex field a as a = \a\e 1 ^ with amplitude \a\ and phase cf>. 
Multiplying Eq. by a*/| a | results in separate differ- 
ential equations for \a\ and <j>: 



a\ — (i>=UNL(a)(- 
dr 



/ Xo+) ~ \/Xo^sm(f) ( 6 ) 
= Vxo~^cos(/>- \a\ (7) 



Eq. J7J) implicates that the amplitude \a\ adjusts expo- 
nentially to y/Xo— cos determined by the instantaneous 
value of </>. This happens at the fastest time-scale avail- 
able inside the cavity, the decay time of the intra-cavity 
field. Hence, any changes of the intra-cavity field on 



slower time-scales are governed by the evolution of <j>. 
This lets us adiabatically eliminate Eq. (J7J by inserting 
\a\= ^/X(7rcos</> into Eq. ©, which leads to 
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Tr 1 



UN 



x/Xo- Xo+ 



L{4>) Xo- cos</>- 



Xo+ 



(8) 



L(<p) 



exp (-£ax ^ . / ; 8 7 ] 
i ^ •suj. cj v Y VXo+Xo- COS0 J 



The localization factor L(<p) is rescaled in terms of 
quantities which remain constant during the time evo- 
lution of the system, i.e. the axial and radial Boltzmann 
factors £,ax and £ ro d, the recoil frequency lur the mode ra- 
dius wo and an effective axial vibrational frequency which 
is a measure of the total power directed to the cavity. 
More precisely, is the axial vibrational frequency cor- 
responding to the optical potential that would arise for 
symmetric pumping with no atoms inside the cavity, i.e. 
if a+ = a_ =£+=£_. 

The dynamical properties of Eq. JSJ) can be analyzed 
in terms of its steady state solutions. The steady state 
values of <f> are obtained as the zeros of the the right hand 
side of Eq. JSJ while the corresponding values of \a\ follow 
directly from Eq. Q . It is particularly instructive to plot 
these steady state values versus the interaction strength 
U N, which can be tuned experimentally. In Fig.|2this is 
shown for different degrees of pumping asymmetry, i.e., 
Xo- = 51%, 50%, 48%, 43% and 38%. In this plot we use 
our experimental data for t; ax , £, ra d, ^v, k and Wq. For 
values of xo- < 43% only one stable solution exists and 
the amplitude and phase of the intra-cavity field a are 
well defined. The situation changes with increasing xo-- 
For higher values, e.g. xo- = 48% a strong bistability 
occurs for an interaction strength of 2 < UN < 3. The 
upper and lower branch with a positive (negative) slope 
for the phase (amplitude) are stable, whereas the mid- 
dle branch is unstable. Experimentally, the interaction 
strength is readily tuned by a change of the particle num- 
ber N. If one increases UN starting from low values the 
system will follow the upper branch of the 48%-trace in 
Fig. [2K) until the turning point at UN w 3 is reached, 
where the intensity suddenly drops to almost zero. On 
the other hand, when we reduce the number of particles 
starting from high values, the intensity abruptly jumps 
from small to high intensity around UN ~ 2. This hys- 
teresis feature is characteristic for bistable systems. If 
we approach symmetric pumping (i.e. xo- approaches 
50% the upper bound of the bistability range (the re- 
gion on the x-axis between the jumps) moves further out 
towards infinity whereas the lower bound only slightly 
increases above two. For xo- > 50% the upper bound of 
the bistability range equals infinity, i.e., the low intensity 
branch cannot be reached by tuning UN. In this case a 
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FIG. 2: Steady state solutions of Equations lRjl7H for the in- 
tensity (a) and the phase (b) of the unlocked intra-cavity 
field for X o- = 51%, 50%, 48%, 43% and 38%. For X o- < 43% 
only one stable solution exists, while for larger values bista- 
bility occurs for UN larger than some value around two. The 
arrows indicate for the 48%-trace, where a jump from one to 
the other stable solution can occur, if the interaction param- 
eter is tuned. 



stable optical lattice can be obtained at any value of the 
interaction strength. 

In our analysis we have neglected non-adiabatic as- 
pects of the atomic motion although this is only justi- 
fied for the axial degrees of freedom which are well con- 
fined. Motion along the radial directions occurs on a far 
slower time scale, and the adiabatic approximation does 
not hold. A typical non-adiabatic reaction of the radial 
motion to sudden changes of the potential depth is the ex- 
citation of breathing oscillations. Such oscillations would 
come along with a small oscillatory change of the effec- 
tive interaction strength. If the system operates deeply 
inside the bistability regime we may immediately infer 
from Fig.|21that small changes of UN with an amplitude 
well below the extension of the bistability range do not 
significantly change the intra-cavity intensity. The most 
drastic reaction of the intensity to a change of UN occurs, 
if the system is close to the frontier between stable and 
bistable operation, as in the 43% trace of Fig. [5] In this 
case we in fact observe oscillations of the intensity which 
are not predicted by our adiabatic model and which are 
subject of Sec. IVII 




FIG. 3: Sketch of the double-MOT system with the resonator 
placed in the vacuum chamber. 



IV. EXPERIMENTAL SETUP 



Preparation of ultracold 85 Rb atoms is established 
with a standard double-MOT (magneto-optical trap) 
setup as depicted in Fig. [3] In a first vacuum cham- 
ber a conventional 3D-MOT collects atoms from a Rb 
dispenser source. One of the retro-reflecting mirrors is 
provided with a hole, such that an atomic beam is ex- 
tracted with a flux of 5 • 10 8 s _1 and a mean velocity of 
lOm/s. This beam loads a second magneto-optical trap 
placed in the main UHV chamber with a pressure below 
2 • 10~ 10 mbar. Here we typically trap 5 • 10 8 atoms at a 
temperature of 100 /iK. The implementation of the MOT 
coils inside the vacuum allows us to switch the magnetic 
field very rapidly. Contamination emerging via heating 
of these coils limits the lifetime of the trap to about 2 s. 

The experimental implementation of the resonator is 
also sketched in Fig. |3| In a triangular setup two curved 
high reflecting mirrors (1.5 ppm transmission) and the 
flat incoupling mirror (23 ppm transmission) form the 
high finesse ring cavity with a round trip path length of 
97mm and a free spectral range of 3.1 GHz. The beam 
waists wq (1/e 2 mode radius) in sagittal and vertical di- 
rections are 134 /im and 129 /im, respectively. The beam 
splitting optics is included in the vacuum to keep the 
path lengths between the splitter and the incoupling mir- 
ror as short as possible. This is favorable because the 
path length difference influences the spatial phase of the 
intra-cavity standing wave. 

For the trapping light we use an extended cavity grat- 
ing stabilized diode laser with variable detuning relative 
to the D2 transition of 85 Rb. An acoustooptical mod- 
ulator (AOM) serves as a fast switch and an intensity 
regulator for the lattice. We analyze the reflected and 
transmitted light with three avalanche photo diodes. 
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FIG. 4: Relative intensity noise power spectrum, a) High 
frequency part of the spectrum up to 1 MHz, relevant for axial 
heating b) Low-frequency part of the spectrum responsible for 
radial heating. The pronounced peak at « 850 Hz originates 
from the turbo pump. 



STABILIZATION, PARAMETRIC HEATING, 
TRAPPING TIME 



In order to stabilize the laser diode emission to the cav- 
ity resonance, we use a frequency modulation technique 
(Pound-Drever-Hall) with a servo bandwidth of several 
MHz. Similarly as in Ref. |2fJ the fast branch of the feed- 
back control directly applies the PDH-error signal to the 
injection current of the diode after passing a loop filter, 
which compensates for frequency-dependent phase shifts 
on the diode laser chip. With the laser locked to the res- 
onator we can measure the photon lifetime l/2j c = 9.3 /is 
corresponding to a finesse of the cavity of F = 1.8 • 10 5 , a 
cavity resonance line width of 17 kHz and a mean scat- 
tering loss per mirror of 3 ppm. 

A sufficiently fast and efficient stabilization is needed 
to prevent intensity fluctuations of the intra-cavity field 
which lead to exponential heating of the trapped atoms 
|2l) . In particular, in the case of a small cavity resonance 
line width tiny frequency deviations are readily converted 
into intensity fluctuations. Unfortunately, the heating 
time scales quadratically with the trap frequency, which 
typically is on the order of a few hundred kHz for the 
axial direction in a standing wave trap pi) : 



W 



;laW + - lr W 7a , r = S(2v*, r ) (9) 



Here W is the mean kinetic energy and v ar are the 
trap frequencies of the atoms in axial and radial direc- 
tion, respectively. S(2v) denotes the one-sided power 
spectrum of the fractional intensity noise. The heating 
rate involves the noise intensity at the second harmonic 
of the vibrational frequency because parametric excita- 
tion is the main source of heating. Here, we assumed that 
the atomic sample is thermalized resulting in the same 
mean kinetic energy for all dimensions. Nevertheless, the 
heating rates for the radial and axial directions may be 
significantly different. We have measured the intensity 
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FIG. 5: a) Temporal evolution of the trap population deter- 
mined by fluorescence imaging. The solid line represents a fit 
to a model which accounts for two-body losses, b) Tempera- 
ture evolution with time. The solid line shows a fit of Ea llUI 
to the experimental data 



fluctuations S{2v) by analyzing the transmission through 
one of the high reflectors. The resulting power spectrum 
of the relative intensity noise is shown in Fig. For typ- 
ical trap parameters f a =350kHz and t^ r =450Hz we find 
S(700kHz)=1.5-10- 13 Hz- 1 and S(900Hz)=3-l(T 9 Hz" 1 
resulting in a heating time (time for an e-fold temper- 
ature increase) 77, = (l/37 + 2/37 r )~ 1 exceeding 20 s. 
This indicates that heating due to intensity noise should 
not be a limiting factor for experiments with trapping 
times up to several seconds. The decrease of the inten- 
sity noise with 20 dB per decade starting at ~17 kHz is 
not surprising, because the resonator acts as a low pass 
filter due to the long photon lifetime. Obviously, the 
small cavity line width only complicates the stabilization 
of the intra-cavity field with respect to low noise frequen- 
cies, while high noise frequencies are strongly suppressed. 

Our noise analysis promises trapping times on the or- 
der of a few seconds limited merely by background colli- 
sions. In order to confirm this prediction we have tuned 
the optical lattice 7nm below the D2-line of 85 Rb, where 
the collective interaction parameter U N is well below 
unity. In this weak coupling regime we have conducted 
life time and temperature measurements. We couple a 
power of 60 fiW into each running wave mode, which is 
enhanced by a factor of about 70000 to produce a trap 
depth of 350 /iK corresponding to axial and radial fre- 
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quencies of 350 kHz and 450 Hz, respectively. 

We typically load 4 • 10 6 atoms into the lattice with 
an initial temperature of about 100 /iK. In Fig. EH) we 
plot the observed trap population versus time. It shows a 
non-exponential two-body decay at the beginning with 
an exponential decay time of 1.7 s representing the main 
contribution after 500 ms. The solid line is a fit according 
to a standard trap decay model including a two-body loss 
term. The two-body losses are explained by evaporation 
arising due to the relatively high density of 9 • 10 11 cm~ 3 
and the large elastic scattering cross section of 85 Rb. The 
expected reduction of the temperature T(t) is shown in 
Fig. |3>). According to a simple theoretical model pre- 
sented in Ref. Q the temporal temperature evolution fol- 
lows 

T(t) = T(0) (i-e^i-e-r*)), 

^dbf 7 - 1 - (10) 

where j3 is the two-body loss coefficient, 7 the expo- 
nential decay due to background losses and po the peak 
density in the potential wells. W denotes the mean ki- 
netic energy removed by an evaporated particle. The 
solid line of Fig. EH>) shows a fit to the experimental data 
by means of Eq. with e = 0.23 used as a fitting pa- 
rameter, while (3 and 7 are taken from Fig. EH) and po is 
measured directly. 

According to Fig. Et>) we do not observe heating at 
least for the first four seconds of the time evolution that 
we can follow before the trapped sample disappears. This 
lets us determine a lower bound for the heating time of 
about 100 s which is a factor of 4 larger than the value of 
Th taken from the noise measurements. This is not sur- 
prising since the temperature is determined for the radial 
direction, whereas the main contribution to heating orig- 
inates from the axial direction. The heating rate along 
the cavity axis is almost two orders of magnitude higher 
as compared to radial heating. While during the first 
500 ms thermalization keeps the axial and radial temper- 
atures locked, this is not so for later times, where the 
density has decreased beyond the collisional regime and 
radial heating should slow down significantly. 

VI. STRONG COUPLING REGIME 

The strong coupling regime is accessed if the collec- 
tive interaction strength Ao N exceeds the decay time of 
the intra-cavity field j c (UN > 1). In order to keep the 
spontaneous scattering time of several ms long enough 
to avoid significant heating, we operate the lattice at 
0.7 nm detuning. This yields a light shift per photon 
A =0.091s _1 . Hence with a few 10 6 atoms trapped in 
the lattice we are able to reach values for UN up to five, 
well within the strong coupling regime. In the following 
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FIG. 6: Scaled intensity \- f° r different values of xo-- The 
MOT beams are shut off at t — ms. 



experiments we have typically applied 5 fiW laser power 
for each travelling wave mode yielding a trap depth of 
about 800 /iK for symmetric pumping and in absence of 
atoms. The corresponding axial vibrational frequency 
of 500 kHz is sufficiently larger than the 17 kHz cavity 
bandwidth, i.e., adiabaticity of the axial motion is a well 
justified assumption. A typical experiment proceeds in 
three steps. We superimpose the MOT upon the optical 
lattice for several seconds, before the MOT light is shut 
off and the atoms remain trapped in the lattice. Finally, 
the lattice is extinguished and the atoms are given some 
time to expand ballistically before a fluorescence image 
of the sample is taken. From the time of flight the ra- 
dial temperature of the sample is determined. Without 
ballistic expansion the fluorescence images give informa- 
tion on the spatial distribution of the trapped atoms. 
We can continuously tune the ratio of powers in the two 
running wave modes from the case of symmetric pump- 
ing xo- = Xo+ to one-sided pumping. We lock the 
laser frequency to only one of the modes, i.e. x+7 an d 
keep its phase and amplitude at a constant value. The 
power leaking out from the unlocked mode through one 
of the high-reflectors is detected and used to determine 
the scaled intensity X-i which corresponds to \a\ 2 in the 
theoretical model in Sec. 11111 

In Fig. El we show the time evolution of the scaled 
intra-cavity intensity X- f° r different types of pumping. 
The MOT is terminated at t = in this figure. For sym- 
metric pumping (xo±=50%) a stable lattice is formed 
with equal intensities in each travelling wave mode irre- 
spective of the number of atoms inside the cavity. For 
values of xo- below 50% the situation changes drasti- 
cally. For an asymmetry of only 1% in favor to the locked 
mode the intra-cavity intensity X- drops to about 10 
percent during the MOT phase and to almost zero af- 
ter shutting off the MOT beams. Subsequently, a slow 
increase is observed before at about 20 ms the intensity 
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FIG. 7: Theoretical simulations of the experimental observa- 
tions of Fig. 
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suddenly jumps back to nearly the value xo- (the value 
expected in absence of atoms) with a rise time below a 
few hundred microseconds. If the power of the unlocked 
mode is further reduced, the time duration before the 
jump occurs increases as well as the rise time until the 
jump completely vanishes for values xo- < 40%. At a 
value of (xo-=43%), where the jump begins to level off, 
strong oscillations of roughly 1 kHz are observed corre- 
sponding to twice the radial vibrational frequency. The 
traces for asymmetric pumping in Fig. are observed for 
Xo- adjusted to 49%, 46%, 43%, 36%, 33%, 26%, and 18% 
respectively. The corresponding values of the initial in- 
teraction strength U N(t = 0) sa 4.48, 4.25, 4.01, 3.54, 
3.30, 2.95, 2.48 are carefully determined by measuring the 
initial particle number N(t = 0) via fluorescence detec- 
tion. The observed decrease of N(t — 0) with decreasing 
Xo- arises because the capture efficiency decreases with 
the lattice well depth. 

In our experiments the value of the interaction strength 
necessarily decreases with time due to trap loss. The 
corresponding decrease of U N in connection with the 
bistability plot of Fig. Ol) explains the observations of 
Fig. As time proceeds in Fig. H3 we move from right 
to left in Fig. [2K) starting on one of the low intensity 
branches in the lower right corner. Hence, depending on 
the value of xo- we encounter a sudden or soft increase 
of intensity, depending on whether we travel on a curve 
with bistable or stable character. We also adjusted values 
of xo- above 50%. In this case X- initially drops to a 
value close to 50% independent of the value of xo- an d 
gradually recovers to xo-- This behavior is understood 
by similar arguments based on Fig. 

We have simulated the time evolution observed in 
Fig. El by means of Eq. The values of |arj| are taken 
directly from the observations of Fig. a) for t = 0. The 
values of r\ ax — 0.5 and r\ ra d = 0.3 are determined by 



FIG. 8: a) Observed intensity of the unlocked mode plotted 
versus time. During the time interval 0ms<t< 2 ms the total 
power coupled to the cavity is reduced by a factor of 2. b) 
shows a numerical simulation based on the adiabatic model. 
The intensity X- (solid line) and the corresponding phase <f> 
(dotted line) are shown. 



temperature measurements with an uncertainty of about 
0.1. The difference in the radial and axial directions origi- 
nates from different trap depths for both directions, since 
the contrast of the interference pattern in axial direction 
depends on the degree of pumping asymmetry. The de- 
crease of N(t) with time is measured and modelled as 
described in Sec. [V] The theoretical simulations shown 
in Fig. [71 reproduce our experimental traces very nicely. 
Not only the general behavior of the jump feature, but 
also the time of the jump to occur, is accurately matched. 
For the initial interaction strength UN we used the values 
2.38, 2.23, 2.15, 1.75 for xo- being 49%, 46%, 43% and 
36%, which fall within a few percent of those determined 
for the corresponding experimental traces, however re- 
duced by a common scaling factor 1.89. The need for 
this factor is not surprising, because in our atom number 
measurements up to a factor two uncertainty should be 
expected for the absolute values, while relative values are 
on the few percent level. 

In the strong coupling regime our system shows an 
anomalous response to a change of the total intensity Jo 
coupled to the cavity while keeping the same pumping 
asymmetry. This is shown in Fig.|8Ji) for xo- = 45%. In 
this plot we have reduced Jo by a factor two at t = ms 
and doubled it again at t = 2 ms. As a first immedi- 
ate reaction to the reduction X- rapidly drops on a time 
scale given by the cavity decay time 7c" 1 , as might be 
expected. However, this is counteracted by an approach 
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FIG. 9: Bistable switching of \- during MOT loading for 

Xo- = 49% 



of an increased steady state value on a slower time scale. 
When the old power level is reestablished, after a tran- 
sient increase, X- drops back to nearly its original value. 
A numerical simulation based on our adiabatic model of 
Sec. reproduces the experimental findings surprisingly 
well except the an unexplained additional ripple at about 
3 kHz, which amounts to about six times the radial vibra- 
tional frequency. The calculations also show that during 
the drop of Iq the phase shift of the unlocked mode with 
respect to the locked mode is reduced. Therefore the ef- 
fective wavelength of the unlocked mode is shifted closer 
to resonance and the intensity increases. 

The bistable character of the atom-cavity system can 
also be observed at work during MOT loading, yielding 
plots as that shown in Fig. 03 where x~ switches between 
two steady state solutions (xo- = 49%). Starting at t=0 
at the high intensity level and hence with a deep lattice 
with comparably low atom number N, the MOT tends 
to increase N until the interaction parameter exceeds the 
critical value and a jump occurs. Now the intensity is 
low and hence the lattice depth, while the temperature 
remains the same. This yields a decrease of the loading 
rate and thus a reduction of N until the system jumps 
back to the previous state. 

In order to control the performance of the locking the 
light reflected from the cavity originating from the locked 
travelling wave mode is monitored together with the light 
of the unlocked mode transmitted through one of the 
high reflectors . This allows to verify, that despite of the 
rapid changes in time observed for X- and c/>, the intra- 
cavity intensity of the counterpropagating locked mode 
remains well behaved. This is illustrated in Fig. for 
Xo-=45%. While in the lower trace rapid time evolution 
is observed the upper trace remains nicely constant. The 
residual structure seen in the upper trace is explained by 
imperfect separation of the contributions from the two 
modes due to limited quality of the polarization optics 
used. 
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FIG. 10: The intra-cavity intensities X- an d X+ °f the un- 
locked mode (lower trace, recorded in transmission) and the 
locked mode (upper trace, obtained from the reflected light) 
are plotted versus time for xo-=45%. 



VII. NON— ADIABATIC MOTION 

The adiabatic model excellently explains the experi- 
mental findings regarding the intra-cavity intensity ex- 
cept for the pronounced oscillation feature observed in 
the xo- = 43%-trace of Fig. EJ As has been discussed 
at the end of Sec. 1 1 1 1 1 the assumption of adiabaticity is 
not well justified for the radial degrees of the atomic mo- 
tion. We came to the conclusion that for fast changes of 
the lattice well depth breathing oscillations should be ex- 
cited which would produce corresponding oscillations in 
the intra-cavity intensity, if the system is operated near 
the frontier between the stable and the bistable regime, 
i.e., in the 43% trace of Fig. [21 In fact this oscillation 
is observed in Fig. [SJ An expanded version is shown in 

Fig. mi). 

As an experimental test of our interpretation in terms 
of radial breathing oscillations we have measured the mo- 
mentum and position spread of the atomic ensemble in 
the radial direction during the observed intensity oscil- 
lations, finding the behavior shown in Fig. Illfal. The 
black rectangles show the radial momentum spread of 
the atomic sample determined by time-of-flight measure- 
ments, while the open circles show the radial spread di- 
rectly observed via in-situ images of the atoms in the 
lattice. The solid and dashed lines are trigonometric fits 
with 7T phase delay, which confirm the expected anti- 
cyclic behavior. In Fig.lllb) we plot the frequency of the 
observed intensity oscillations versus the potential well 
depth confirming the expected square-root dependence. 

In order to include non-adiabatic aspects in our theo- 
retical description we use the full set of 6N+2 equations 
of motion of Eq. J2J and Eq. (i.e., 6N equations for 
the atomic positions and momenta and two equations for 
the amplitude and the phase of the unlocked intra-cavity 
mode). Since these equations are coupled and non-linear, 



10 




500 1000 1500 

Time [us] 




Well Depth [jxK] 

FIG. 11: (a) Oscillation of the atomic position spread (cir- 
cles) and momentum spread (rectangles). The solid and 
dashed lines are trigonometric fits with tt phase delay, (b) 
Oscillation frequency plotted versus the well depth. The solid 
line shows the expected square root dependence. 
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FIG. 12: a) Experimental observation of radial breathing 
oscillations, (b) Simulation of the oscillations of (a) by means 
of solving the complete set of equations of motion for hundred 
atoms with an upscaled interaction strength. 



the simulation of all 10 6 atoms is beyond our computa- 
tional capacities. Therefore, in our calculations we re- 
duce the number of atoms to one hundred and work with 
an increased light shift per photon, such that the interac- 
tion strength acquires values which compare to the exper- 
iments. The artificially increased light shift per photon 
comes along with a correspondingly increased light shift 
acting on each atom. Hence, in order to maintain the 
potential well depth at the level used in the experiments, 
we work with correspondingly decreased incoupled inten- 
sities. 

The full model for one hundred atoms very accurately 
reproduces our calculations based on the adiabatic model 



shown in Fig. [7| except for the fact that in the vicinity 
of Xo- ~ 43% additional oscillations arise. This confirms 
once again that the assumption of adiabaticity of Sec. 
is well justified for the axial degrees of freedom. We can 
quantitatively reproduce the frequency of the squeezing 
oscillation as is shown in Fig. ^]b). The experimental 
data are taken for UN w 3.5 and an axial vibrational 
frequency of 500 kHz, whereas for the theoretical curve 
UN — 2.0 and 550 kHz was used. Applying the scaling 
factor 1.89 compensating for our systematic overestima- 
tion of the atom number similarly as in Sec. IVII the mea- 
sured interaction strength coincides with the theoretical 
value within 10%. The discrepancy in the vibrational 
frequencies lies within our observation uncertainty of ca. 
10%. The non-adiabatic simulations of the system also 
reproduce the small frequency decrease in time observed 
in Fig. [Hi). 



VIII. CONCLUSIONS 



We conclude that optical lattices formed inside high fi- 
nesse cavities open up an interesting new regime charac- 
terized by collective interactions significantly contribut- 
ing to the atomic dynamics. This regime can be experi- 
mentally realized even far from an atomic resonance such 
that the interaction is entirely of dispersive nature. Al- 
though the combination of a high finesse with a large 
mode volume yields a small cavity resonance bandwidth, 
the experimental challenge of preventing intra-cavity in- 
tensity fluctuations can be handled. In the ring cav- 
ity studied here, specific parameter ranges are identified 
which allow to operate a stable lattice independent of 
the number of trapped atoms. Other regimes are found 
where dispersive optical bistability accompanied by self- 
induced breathing oscillations occur. The spatial phase 
of the lattice is not pinned by the phases of the incou- 
pled laser beams, but rather determined by the strength 
of the atom-cavity interaction. The bistable behavior 
arises for asymmetric pumping and can be understood in 
terms of an adiabatic approximation for the atomic mo- 
tion, whereas the general equations of motion must be 
considered to model the observed breathing oscillations. 
In this article we have only discussed selected aspects 
of the system dynamics. Other interesting phenomena 
could be studied, as for example collective atomic recoil 
lasing (CARL) |22j, which has been recently observed for 
unidirectional pumping of the ring cavity Q ■ The strong 
coupling regime might also be utilized for implementing 
novel laser cooling schemes, which rely on cavity-tailored 
coherent scattering, for example as described in Ref . ^3 • 
Such schemes are highly desirable since they promise to 
extend laser cooling to new species and to operate in a 
density regime not yet accessible. 
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